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ABSTRACT 

This paper presents simple analytic approximations to the linear power spectra, 
linear growth rates, and rms mass fluctuations for both components in a family of 
cold+hot dark matter (CDM+HDM) models that are of current cosmological interest. 
The formulas are valid for a wide range of wavenumber, neutrino fraction, redshift, and 
Hubble constant: k^lOh Mpc"\ 0.05 ^ ^ 0.3, < z ^ 15, and 0.5 ^ h ^ 0.8. A 
new, redshift-dependent shape parameter = a^^^ili^h"^ is introduced to simplify the 
multi-dimensional parameter space and to characterize the effect of massive neutrinos 
on the power spectrum. The physical origin of lies in the neutrino free-streaming 
process, and the analytic approximations can be simplified to depend only on this 
variable and Linear calculations with these power spectra as input are performed 
to compare the predictions of ^ 0.3 models with observational constraints from the 
reconstructed linear power spectrum and cluster abundance. The usual assumption of 
an exact scale-invariant primordial power spectrum is relaxed to allow a spectral index 
of 0.8 ^ n < 1. It is found that a slight tilt of n = 0.9 (no tensor mode) or n = 0.95 
(with tensor mode) in ~ 0.1 — 0.2 CDM-I-HDM models gives a power spectrum 
similar to that of an open CDM model with a shape parameter P = 0.25, providing 
good agreement with the power spectrum reconstructed by Peacock and Dodds (1994) 
and the observed cluster abundance. 



Subject headings: cosmology : theory - dark matter - elementary particles - large-scale 
structure of universe 



1. Introduction 

The linear power spectrum of a given cosmological model serves as the basic input for most 
linear calculations of large-scale structure and numerical simulations of gravitational collapse. 
The computation of the power spectrum involves numerical integration of the linearized, coupled 
Einstein, Boltzmann and fluid equations that govern the evolution of the metric perturbations and 
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the density fields of different particle species. The calculations become increasingly complicated as 

more particle species are included, with those for massive neutrinos being the most time-consuming 
due to their time-dependent energy-momentum relation and nonzero anisotropic stress (Ma &; 
Bertschinger 1995 and references therein). 

The Boltzmann integration is necessary for obtaining highly accurate power spectra, but 
when an accuracy of ^ 10% is all that is desired, analytic approximations can be very useful 
and illuminating. The power spectra for O = 1 CDM models (with or without a cosmological 
constant) have been computed and various fitting functions have been given (e.g., Bardeen et al. 
1986 for zero-baryon models; Holtzman 1989; Efstathiou, Bond, & White 1992). But to date, 
no complete spectra have been published for both components in the CDM-I-HDM models for 
all cosmologically interesting ranges of parameters. Specific models have been examined, with 
various approximations made in the Boltzmann calculation. Holtzman (1989) compiled fitting 
formulas for the present-day baryon transfer function for 0,,^ = 0.1 and 0.3. Van Dalen &; Schaefer 
(1992) tabulated fits to the previously computed transfer functions (Schaefer, Shafi, & Stecker 
1989) for the density-weighted power spectrum at z = 10 for discrete values of between 
and 0.53. Klypin et al. (1993) studied the ^li, = 0.3, Qh = 0.1 model. Pogosyan & Starobinsky 
(1995) provided a fitting formula for the linear gravitational potential. Although the total power 
spectrum is useful for many linear calculations, the separate CDM and HDM spectra and the 
growth rate are essential, for example, for generation of initial conditions in numerical simulations. 
Moreover, all the Boltzmann calculations above either assumed zero baryons, which can lead 
to discrepancies as large as 25% for models with 5% baryons (see §2.2 below), or ignored the 
anisotropic stress and the higher moments in the distribution function of the massive neutrinos. 

Section 2 of this paper presents analytic approximations to the linear growth rate of the 
density field, the CDM and HDM power spectra, the density-weighted power spectrum, and 
the rms mass fluctuations in spatially flat 0,^, ^ 0.3 CDM-I-HDM models. The approximations 
generally differ by less than 10% from the numerically integrated results for the ranges k ^ lOh 
Mpc^^, < z ^ 15, and 0.5 ^ h ^ 0.8. The power spectra used in the fit are computed from 
the Boltzmann code discussed in Ma &; Bertschinger (1995), which includes a careful treatment 
of the massive neutrino phase space and uses 50 Legendre-modes in the angular expansion of the 
distribution function. The calculations include CDM, baryons, photons, and one massive and 
two massless neutrino species. The range < 0,^ ^ 0.3 is chosen since it spans the standard 
(albeit troubled) pure CDM model {Qi, = 0) and the previously preferred Oj, = 0.3 model in the 
literature. Late galaxy formation due to excessive neutrino free-streaming invalidates the fij, = 0.3 
model (Mo & Miralda-Escude 1994; Kauffmann & Chariot 1994; Ma & Bertschinger 1994; Klypin 
et al. 1995), leaving < 0,,^ < 0.3 the only open window for the standard CDM-I-HDM models. 
Numerical simulations of COBE- normalized CDM-I-HDM models indicate that for an exact 
scale-invariant power spectrum (n = 1), a value of 0,^, ~ 0.2 is needed to simultaneously satisfy the 
high-redshift constraint from damped Lyman-a absorbers and the zero-redshift constraint from 
galaxy pairwise velocity dispersions (Ma &; Bertschinger 1994; Ma 1995). 
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In §3, linear theory, with the power spectra above as an input, is used to explore the Qrms-PS! 
n, h, and parameter space. The task of gathering observational constraints and finding a 
concordance region in various slices of this multi-dimensional parameter space has been performed 
by van Dalen & Schaefer (1992), Liddle & Lyth (1993), Pogosyan & Starobinsky (1993, 1995), 
and Liddle et al. (1995). The emphasis here is on physical understanding and explicit illustration 
of the parameter-dependence of the growth rates, the power spectra, and cluster abundance. The 
effect of a small tilt (n ^ 0.8) in the primordial power spectrum on the cluster abundance is also 
examined. Setting the spectral index n to unity in the standard CDM and CDM-I-HDM models is, 
after all, only a simplification. Most inflationary models predict nearly scale-invariant spectra for 
the density fluctuations that have an additional logarithmic dependence on the wavenumber k. The 
primordial power spectrum is therefore never exactly scale-invariant, although a /c-independent, 
n ^ 1 spectral index approximates the logarithmic dependence well (Crittenden & Steinhardt 
1992). The intention here is to explore small deviations from the canonical n = 1, but not to study 
the models with large tilts that would require finer tuning of the inflaton potentials. 



2. Analytic Approximations 

The presence of massive neutrinos in CDM-I-HDM models complicates the computation of the 
power spectrum because the initially relativistic neutrinos become non-relativistic at ksTi, ~ rriuC^- 
This is reflected in the time-dependent energy- momentum relation: = q^c^ + a'^mlc^, where e 
is the comoving energy, a is the expansion factor, and q denotes the magnitude of the conjugate 
momentum, which is a constant for free-streaming neutrinos in the absence of metric perturbations. 
As a result, although the g-dependence in the phase space distribution can be integrated out for 
massless particles in the Boltzmann calculation, the distribution of massive neutrinos must be 
computed on an additional grid for the momentum space. This lengthens the numerical integration 
time considerably, especially for large wavenumbers (see, e.g.. Ma & Bertschinger 1995). 

The presence of massive neutrinos also considerably complicates the procedure for fitting the 
power spectrum because the free streaming of the neutrinos retards the growth of the density 
fields, but only below the free-streaming scale. The amount of suppression depends on the redshift, 
the length scale, the Hubble constant, and the neutrino mass (or equivalently, ^u)- As a result, 
the linear power spectrum does not simply grow as on all length scales as in the standard 
CDM model, and the analytic approximations will depend on Q^y, h, and a in addition to the 
wavenumber k. The simple fitting procedure to a function Pc{k) of a single variable in the CDM 
model therefore becomes a multi- variable problem when massive neutrinos are present. 

Figure 1 shows the numerical results from our Boltzmann code (available at 



http: / / arcturus.mit.edu / cosmics / ; see Bertschinger 1995 and Ma & Bertschinger 1995 for 
descriptions) for the present-day (a = 1) power spectra for the pure CDM and four CDM-I-HDM 
models. The spectral index is taken to be n = 1, and Hq = 50 km s^^ Mpc^^. One would perhaps 
first adopt the brute force approach of fitting the cold and hot spectra Pc and Pu directly for the 
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ranges o{ k,Ui,, a, and Hq of interest. However, Figure 1 shows that unlike Pc, Pi, does not vary 
systematically with Q,y for k ^ 0.5h Mpc"-*^. This behavior makes it difficult to envision a simple 
functional form to approximate P^,. 

One solution is to fit the ratios of the power spectra rather than the spectra directly. Figures 2 
and 3 show the behavior of three functions /, g, and H that are defined by 



The function f{k, a, represents the growth rate of the CDM density field, g{k, a, $7^) represents 
the ratio of the CDM power spectrum P^ in CDM+HDM models to that in the pure CDM model, 
and 7i{k, a, ^lu) represents the ratio of the power in the hot and cold components in a given model. 
The merit of using these functions is that they all vary systematically with fl^, k, and a, as shown 
in Figures 2 and 3. It is therefore much easier to formulate an elegant ansatz for the analytic 
approximations and obtain greater accuracies for the fits. 

There are still four parameters k,U,y, a, and Hq to be considered. Further simplification 
is achieved by identifying physical processes that introduce special length scales into the power 
spectrum. It is useful to first recall that the power spectrum for a pure CDM model exhibits a 
break at keq, the horizon size at matter-radiation equality. Since keq cx flmh'^ for a model with 
a total matter density ilm, and the observable wavenumbers are in units of h Mpc^^, the shape 
parameter F = Q^h can be used to characterize this single feature in the pure CDM power 
spectrum (Efstathiou et al. 1992). The same break should also exist in the power spectra of 
CDM+HDM models since the density perturbations in both components cannot grow if they enter 
the horizon in the radiation-dominated era. In this paper, = 1, so F = /t. 

For the CDM+HDM models, the neutrino free-streaming process introduces a second length 
scale in the power spectrum. The scaling of the free-streaming distance can be understood in terms 
of the (comoving) neutrino Jeans wavenumber kf^ = AirGpa? /v'^^^, where Umed is taken to be the 
median neutrino speed in the Fcrmi-Dirac distribution. (See Bond & Szalay 1983 for a discussion 
of Jeans lengths for collisionless particles.) For k < the density perturbation in the neutrinos 
is Jeans unstable and grows unimpeded in the matter-dominated era; for k > fcfg, the density 
perturbation decays due to neutrino phase mixing. When the neutrinos are relativistic, Umed ~ c , 
and the free-streaming distance is approximately the particle horizon, which scales as k{s{a) oc 
(in the radiation-dominated era). The neutrinos become non-rclativistic at a ~ ?>kBTQy/my(? 
(where Tqi^ is the neutrino temperature today), after which Wmcd oc kBT^y/arriyC. The relations 
777,1/ oc and p oc a~^h? then imply fcfs(a) oc a}/^Q.i,h^ . The decrease with time in the 

free-streaming distance of non-relativistic neutrinos allows neutrinos to fall into the CDM potential 
well and to grow with the CDM on increasingly smaller scales. This effect can be seen by 
comparing / and H in Figure 2 and 3, where both functions remain at unity to a larger A; at a = 1 



f(k,a,n„) 



1 dlogPc{k,a,Q„) 



Pc{k,a,Q.y) 
P„{k,a,Q,y) 



2 d\oga 

g{k, a, 0,1,) Pc{k, a, 0,^, = 0) 
TC{k, a, rij,) Pc{k, a, 0„) . 
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than a = 0.1. 

Based on the reasoning above, one can simphfy the dependence on the parameters a, and 
h by introducing a second shape parameter, 

= aV2o,/i2 , (2) 

that characterizes the effect of decreasing free-streaming distances on the power spectrum. This 
reahzation was crucial for the simphfication of the functional forms and the improvement of the 
fits described below. 

Below I describe the analytic approximations to the linear growth rate, the separate CDM 
and HDM power spectra, the density-weighted power spectrum, and the rms mass fluctuations. A 
summary of the equations and the fitting coefficients is given in Table 1. 



2.1. Growth Rates 

In the standard CDM model, the CDM density field 5c grows as the expansion factor on all 
scales; therefore / = d\og5 / dloga = 1. The massive neutrinos in CDM+HDM models introduce 
an additional length scale, the free-streaming distance, below which the density fluctuations are 
washed out and the growth rate is retarded. Consequently, the growth rates in CDM+HDM 
models are functions of the wavenumber k, Q.^, and time. This dependence is illustrated in 
Figures 2a and 3a. 

Before searching for an analytical approximation, it is useful to ascertain first the asymptotic 
behavior of the function. Since HDM behaves like CDM above the free-streaming distance, / ^ 1 
as /c — > 0. At large k where the HDM density field is negligible compared to CDM, the time 
evolution of the latter is governed by the linearized fluid equation 5c + a5c/a = l.bH'^a?Q-c5ci 
where the dots denote differentiation with respect to the conformal time r. Since Ha = 2/t in 
the matter-dominated era, the growing solution is easily shown to give (Bond, Efstathiou, & Silk 
1980) 

/oo = f{k ^ oo) = + \ = yi-^n, - i . (3) 

It is interesting to note that 

foo « (4) 

is an excellent approximation to equation (|3|), especially for the range Qc ^ 0.7 studied in this 
paper. It is sometimes convenient to express all dependence in terms of 17,^ = 1 — Qc, and 

/oo « n'c' « 1 - o.Gsnl-^' (5) 

is accurate to better than 99.7% for i^^, ^ 0.3 (for both expressions). 

It is not a coincidence that the form above is identical to the widely-used formula / = for 
a mildly open model with matter density equation (^ is exact for the growth rates in open 
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models if f^c is replaced with and a CDM+HDM model in fact evolves like an open universe 
with an effective density f^c < 1 on scales much below the neutrino free-streaming distance. 
Nevertheless, it should be remembered that while / is scale-independent for open models, it is 
scale-dependent for CDM+HDM models. 



An excellent approximation to the growth rate is given by 



1 + ^7g-6Qi x"^ _ k 

f(x,SZi/) = , X = —— , (6) 

^ ' ^ 1 + ai x«'2 ' T^h' ^ ^ 

where Ti, = a}-/'^Vlj^h? , Q.^ + = 1, and the best-fit parameters are oi = 0.1161 and 02 = 1.363 
for k in units of Mpc~^. The fractional error of the fit relative to the numerically computed values 
is smaller than 0.5% for k ^20h Mpc'^, 0.05 ^ 0^ ^ 0.3, 0.5 ^ h ^ 0.8, and < z ^ 15. 
The functional form is chosen to approach the asymptotic values discussed above. As promised, 
equation (^) depends only on the variable x that characterizes the neutrino free-streaming scale, 
and i^i, (or J7c) via foo- The solid and dashed curves in Figures 2a and 3a compare the exact and 
fitted growth rates at o = 1 and 0.1 in different CDM+HDM models with h = 0.5. The top panel 
of Figure 4 illustrates the perfect scaling with the Hubble constant through k/T^h oc k/h? . The 
shape parameter F associated with matter-radiation equality discussed earlier does not affect / in 
the matter-dominated era of concern here. 



2.2. CDM Povi^er Spectrum 

Since the strategy here is to use the ratios of the spectra, it is essential to have an accurate 
approximation for the power spectrum Pc{k.,a.,Q-y = 0) in the pure CDM model as a starting 
point. I will adopt the functional form of Bardeen et al. (1986) 



ln(l + aiq) 



k 



aiq J [l + a2q + {a3qy + {aiqf + {a5q)^V/^' ^ F/i ' 

(7) 

where F = as discussed above. However, the coefficients ai = 2.34, 02 = 3.89, 03 = 16.1, 04 = 
5.46, and = 6.71 from Bardeen et al. (1986) are accurate for the zero-baryon CDM model only. 
Increasing the baryonic fraction fib decreases the power since the baryonic component cannot fall 
into the CDM potential wells until after recombination. Compared to the standard CDM model 
with = 0.95 and Qj, = 0.05, equation (^) overpredicts the power by as much as 25% at k ^ 0.1 
Mpc"^ (for fixed long- wavelength normalization A). A high-accuracy fit for the Qb = 0.05 CDM 
model can be achieved by modifying the coefficients to ai = 2.205, 02 = 4.05, 03 = 18.3, 04 = 8.725, 
and = 8.0. The fractional error relative to the direct numerical result is smaller than 1% 
for k < 40h Mpc~^. Alternatively, replacing F = Om/i with F = exp{—2^lh)^mh and using the 
Bardeen et al. coefficients in equation (0) can approximate the effect of the baryons (Peacock 
& Dodds 1994). This enables one to explore a wide range of 0,b, but the error is larger for the 
standard i^b = 0.05 model: ^ 10% compared to ^ 1% provided by the new coefficients above. 
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If the temperature fluctuations arise purely from the Sachs- Wolfe effect (Sachs &: Wolfe 1967), 
the normalization factor A in equation (|^) is related to Qrms-PS by 



A = k'-''—[^\ 211^^ = 2.689 X 10^ /i-*Mpc% (8) 

° 207r \HoJ ^ 1^^2.726/ ^ ^ 

where Qis = Qrms-Ps/18 /xK , r2.726 = ro,^/2.726i^ , and ko ^ IHq/2c « 1/600 /i Mpc^^ for 
COBE. There are other corrections to this formula, but it is good to ^ 1% for = \ CDM models 
with /i = 1 and ^ 4% for h as low as 0.3 (Bunn, Scott, & White 1995; White & Bunn 1995). The 
best COBE value for Qrms-PS is also model-dependent and is discussed in §3.1. 

The next step is to relate the CDM power spectra Pc in models with different Vty. Figures 2b 
and 3b illustrate the effect on when CDM is partially replaced by HDM in a spatially-flat model. 
Since the massive neutrinos do not cluster until the neutrino temperature has dropped below their 
mass scale, the CDM component evolves more slowly early on with an effective Q, = \ — Qy < 1 &t 
k above the free-streaming scale . Consequently, the ratio g = Pc{k,Qy)/ Pc{k,Qi, = 0) becomes 
smaller for larger ^1^. 

The large-A; behavior of g can be derived analytically from the definition of the growth rate in 
equation (||). It is g{k oo) oc a^*^-^°°~^\ where foo is given by equation @. Using 1 — /oo oc il.}, '^^ 
as given by equation (^), I find the following functional form to give a good approximation: 

(Ql05 
: : — 6T , ^ = ' xo = x{a = l), (9 

where = a!" I 0.^1? , and the best-fit parameters are 6i = 0.01647,62 = 2.803 x 10^^ 63 = 10.90, 
and 64 = 3.259 for k in units of Mpc^^. It is important to note that x is defined to be kjVy 
instead of kjVyh as in equation (^. This results from the fact that the two shape parameters 
Y = h = 1 in this paper) and = a^/^O^/i^ discussed earlier depend on different powers 
of the Hubble constant, and both parameters affect the shape of 5, which involves Pc from two 
different models. Since the free-streaming wavenumber /cfg is much larger than fcgq at the redshifts 
of concern here {z ^ 20), F oc /i tends to control the scaling with h in the region where g starts 
to deviate from unity in Figures 2 and 3, while F^, oc controls the scaling with h at larger k. 
This mixed dependence on h causes potential complication for the fitting formula. For the range 
0.5 ^ h ^ 0.8 and k ^ Wh Mpc~^ of interest here, however, I find that the scaling x oc k/h'^ 
works fairly well, as illustrated by the middle panel of Figure 4. Therefore, x is modified to 
depend on one less power of h than for the function / to account for the effect of F. The error 
in Pc in equation (^) relative to the numerically computed values is <^ 10% for k <^ Wh Mpc^^, 
0.05 ^n^^ 0.2, 0.5 ^ /i ^ 0.8, and < z ^ 15. For ^ 0.3, the error is stih ^ 10% for h = 0.5 
and ^ lO/i Mpc^^, but it increases to ~ 20% for h = 0.8 at high k. 
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2.3. HDM Power Spectrum 

The HDM to CDM ratio 7{{k,0,u) decreases monotonically with decreasing ili, as shown by 
Figures 2c and 3c. This trend is caused by the longer neutrino free-streaming lengths for smaller 
neutrino masses (hence smaller il-u), which makes Ti. deviate from unity at smaller k. As discussed 
earlier, the dependence of Pu/Pc on k, a, and Hq can be combined in a single variable x that 
characterizes the neutrino free streaming distance. A good analytic approximation is given by 

Pe(A:,a,0,) " "^^^^ " l + C2xl/2+c3X + C4x3/2 + C5x2 ' " T./l ' 

where = a^/^VL^h^, and ci = 0.0015, ca = -0.1207, cg = 0.1015, C4 = -0.01618 and 

C5 = 0.001711. The resulting Pi^ has an error ^ 15% relative to the numerical results for k ^IQh 

Mpc-\ 0.05 <,Vly<, 0.3, 0.5 ^ /i ^ 0.8, and < z ^ 15. 

The bottom panel of Figure 4 illustrates the perfect scaling of 7i with x oc k/h?. The 
complication due to different /i-scalings in the function g (see §2.2) does not arise here because 
the physical process associated with shape parameter V = h introduces breaks in Pc and P^ at the 
same length scale (i.e., at ^eq), and the effect interestingly cancels out in the ratio Ti. 



2.4. Density- Weighted Power Spectrum 

Although the clustering of the cold and hot components is described by the separate power 
spectra Pc and Pu discussed above, it sometimes is useful to have an analytic approximation to 
the density-weighted power spectrum P{k) = {QyPl^'^ + {1 — Qy)Pl^'^}'^ that measures the total 
gravitational fluctuations contributed by the separate components. Here I have assumed that 
CDM and baryons have the same power (i.e., Pc = Pb), which is a very good approximation for the 
range of redshifts and $7^ studied in this paper. The functional form used for the CDM spectrum 
Pc in equation (|9|) works well here, and a good approximation is given by 

CQl-05 
I -\- di x'^^l'^ + x'^^ \ " k 
> ^ = T^' XQ = x{a = \), (11) 
1 + Xq^ / i 1/ 

where Pc{k,a,Cli, = 0) for the pure CDM model is given by equation (|^), and 
di = 0.004321, (i2 = 2.217 x 10"^, ^3 = 11.63, and d^ = 3.317. The error in the resulting 
P relative to the numerically computed values is <^ 10% for the ranges k <^ 20h Mpc~^, 
0.05 ^n^^ 0.3, 0.5 ^ /i ^ 0.8, and < z ^ 15. 

Pogosyan & Starobinsky (1995) introduced a more complicated 8-parameter fitting formula 
for the gravitational potential that gives results similar to ours. The analytic approximation 
provided by Efstathiou et al. (1992), however, does not fit the density- weighted CDM+HDM 
P{k) well. Their formula P{k) = Bk/{1 + [ak + (bk)^/^ + (ckf]"}'^/'' , with a = 6.4/(r/i) Mpc, 
b = 3.0/(r/i) Mpc, c = 1.7/(r/i) Mpc, u = 1.13, and F = 0.2(J7j,/0.3)-°-5 does not take into 
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account the redshift dependence of the shape and does not have the correct dependence on Cli, at 
large k. It produces too Uttle power at 0.01 ^ k <^ 1 Mpc~^ by as much as a factor of ~ 2. On the 
other hand, this formula with T = flmh = 0.5 works fairly well for the fij, = 0.05, h = 0.5 flat CDM 
model. The maximal discrepancy is ~ 18%, compared to 25% in the original Bardeen et al. fits. 



2.5. Mass Fluctuations a 

The linear rms mass fluctuations in spheres of radius R are related to the power spectrum P 

by 

a^iR,n^,a)= —4TTk^P{k,n^,a)W^{kR), (12) 
Jo k 

where W{x) = 3(sinx — xcosx)/a;^ is the tophat window function, and the mass enclosed in 
the sphere is M = 47rpoR^/3, with po denoting the background mass density of the universe. 
Although different filter windows have been used in the literature, and there is no a priori reason 
to prefer a particular window function, Lacey &; Cole (1994) find that the tophat gives the best 
fit to A?^-body simulations with initial power-law P{k). More importantly, the best-fit Sc for a 
tophat is insensitive to the spectral index n. In contrast, they find the best fit Sc to depend more 
strongly on n when the Gaussian window is used. Ma & Bertschinger (1994) also find the shapes 
of the Press-Schechter mass functions computed with the tophat filter to give a better match in 
CDM+HDM simulations. Since a family of models with different spectral shapes at large k is 
being studied here, the tophat window will be the most robust filter. 

Figure 5 shows the present-day a{R, for n = 1 and h = 0.5 models with different Q,^. At 
large R, all curves converge because the same 4-year COBE normalization (Qrms-PS = 18 /J-K) is 
used, and a oc since P{k) oc k" at small k. The small-i? behavior simply refiects the 

effect of neutrino free streaming on the power spectrum. I will again fit to a{Qi, = 0) for the 
pure CDM model and then the ratios a{i^i,)/o'i'^u = 0). An excellent analytic approximation to 
(7(rii/ = 0) is given by 

a{R,n, = 0) = Qish^ — ) -— — 2' x = Rh\ (13) 

where Qis = Qrms-Ps/lS^/C, ai = 0.01359, 02 = 0.05541, 03 = 0.001702, 04 = 0.8032, Rq = 1000 
Mpc, and R is in units of Mpc. The exponent (1 — n)/2 is chosen so that a oc at large 

R. The fractional error is ^ 5% for R = 0.07 - 350 h''^ Mpc, or M = 4 x 10^ - 5 x 10^^ h'^MQ. 
For CDM-I-HDM models, one can use 

a(R,n^) = a(R,n^ = 0)— , x = RT^, , (14) 

1 -I- ai a;"2 ' 

where = ^yh?, 0.^ = 1 — O,^,, and the best-fit parameters are ai = 0.7396, 02 = —0.8927, and 
aa = 5.106. The fractional error is ^ 10%. The result of the fitting is shown in Figure 5 (as dashed 
curves) . 
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As emphasized earlier, the shape of P{k, a, i}^) has a non-trivial dependence on the redshift due 
to neutrino free streaming. However, the mass fluctuations obey a{R,^y,a) = aa{R,Vl^,a = 1) 
to a good approximation, with <10% error for M ^ IO^^Mq (see Figure 5). I therefore will not 
give a time-dependent approximation. If a is required to higher precision, or at a lower mass, one 
can carry out the simple integral in equation (|l^) using the analytic approximation for P{k^ a, Vly) 
in equation ([lT|) . 

3. Observational Constraints 

3.1. Reconstructed P{k) 

The theoretical linear power spectra computed above can be put to the test against the 
observed power spectra of galaxies and clusters, if the effects of non-linear evolution, redshift-space 
distortions, and bias are corrected for. This is not a straightforward task, especially on highly 
nonlinear scales. By limiting the study to the quasi-linear regime [k <, 0.5 h Mpc~^), Peacock & 
Dodds (1994) attempted to reconstruct the linear mass power spectrum from eight independent 
sets of data from optical, radio, and IRAS galaxies and optical and radio Abell clusters. They were 
able to obtain good internal agreement for the estimated power spectrum given specific conditions: 
a significant redshift-space distortion; a scale-independent bias for a given class of objects, with 
relative bias factors for Abell clusters, radio galaxies, optical galaxies, and IRAS galaxies required 
to be 4.5 : 1.9 : 1.3 : 1 to within 6% rms. Among the pure CDM models parameterized by the 
shape parameter T = (see eq. ||7)), the best fit was found to be given by F 0.25. 

Figure 6 compares the final estimate for the linear power spectrum from Peacock & Dodds 
(1994) with those for the standard CDM and three CDM-I-HDM models with h = 0.5 computed in 
the previous section. (The density-averag ed P = {nuPu^^ + {l-^u)Pc^^}^ is used for CDM+HDM 
models.) Three primordial spectral indices n = 1, 0.95, and 0.9 are shown. The bottom two 
panels illustrate the effect of tensor contributions in the tilted models. Whether a model produces 
a non-zero tensor fluctuation depends on the shape of the inflaton potential at the time of the 
first horizon crossing (60 or so e-folds before the end of inflation). Inflationary models, such as the 
extended or chaotic models, generally give a tensor-to-scalar ratio of T/S ~ 7(1 — n) (Davis et al. 
1992), although the so-called natural inflationary model predicts a negligible tensor contribution 
(Adams et al. 1993). Since COBE measures the combined scalar and tensor anisotropics, the 
effect of a non-zero tensor mode, when the tilt is small, is to lower the normalization for the scalar 
power spectrum hy S/{S + T) for a given Qrms-PS- 

An additional effect on the normalization that must be considered is the dependence of 
Qrms-PS on parameters such as n, Qf, and h assumed in a given model. The dependence arises 
because there exist small corrections to the Sachs- Wolfe effect (Sachs & Wolfe 1967) that is the 
main contributor to the angular power spectrum for the temperature fluctuations in the cosmic 
microwave background on the COBE scale. A detailed analysis of such dependence has been 
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carried out by Bunn et al. (1995), who obtain Qrms-PsiiT-) = 21.1 exp[0. 69(1 — n)] /iK from the 
2-year COBE data (for a pure Sachs- Wolfe spectrum). A maximum Hkehhood analysis of the 
same data by Gorski et al. (1994) finds Qj-^^— psin) = (39 + 2n) exp(-0.73n) with quadrupole, 
and Qrms-Ps(^) = (40 + 2n) exp(— 0.74n) without quadrupole. Results from the 4-year COBE 
data have been announced recently (Gorski et al. 1996; Bennett et al. 1996). I therefore will 
use the latest value Qrms-PS = 18 fJ-K for n = 1 (Gorski et al. 1996), and Qrms-PS = 19.2 fiK 
and 20.5 fxK for n = 0.9 and 0.8 (Gorski, private communication), respectively. (Note that 
using (5rms-ps(l) = 18 instead of 21.1 in the Bunn et al. formula gives very similar results: 
Qrms-PS = 19.3^7^ for n = 0.9 and 20.7 /iK for n = 0.8.) The dependence on n for n = 0.8 — 1 
is weak, but it is taken into account in Figure 6. Varying fi^ from 0.01 to 0.1 in the CDM model 
changes Qrms-PS by only about 1% (Bunn et al. 1995), so it will not be considered here. 

As Figure 6 illustrates, the 4- year COBE result is consistent with the amplitude of 
reconstructed power spectrum at small k, but all n = 1 models with < il.,^ < 0.3 normalized to 
this value have too much power at 0.05 ^ k <, 0.2/iMpc^^. A slight tilt of n = 0.9 — 0.95 brings 
the CDM+HDM models into excellent agreement with the data. For a given COBE normalization, 
the power spectrum of a tilted model with a negligible T/S (bottom panel of Fig. 6) has a higher 
amplitude than that of a model with T/S ~ 7(1 — n) (middle panel); therefore a slightly smaller 
n is generally needed to match the reconstructed power spectrum. Figure 6 shows that the four 
highest k points in Peacock & Dodds (1994) favor Q^, ^ 0.1 — 0.2, although Liddle et al. (1995) 
recently questioned the validity of these high-A; points. It is also interesting to note that the 
spectra for the = 0.1 model with n = 0.95 (with tensor) and n = 0.9 (without tensor) are very 
similar to the open CDM spectrum with F = 0.25. 

Figure 7 shows the effect of increasing h to 0.65. The discrepancy worsens for n = 1, but a 
larger tilt of n ~ 0.9 (with tensor) or n ~ 0.8 (without tensor) can compensate for the increased 
power and bring the models with fi^^ ~ 0.1 — 0.3 into good agreement with the data points. 



3.2. as and Cluster Abundance 

The rms linear mass fluctuation on scales of 8h^^ Mpc, as, has played a very special role 
in the literature. Because the observed rms galaxy count on this scale is about unity (Davis & 
Peebles 1983), as was used to define the "bias factor", b = a^^ , between mass and galaxies. It was 
also a common practice to use as for the normalization of the linear matter power spectrum, a 
value undetermined by most theories. After COBE, the rms quadrupole inferred from the cosmic 
microwave background anisotropy is often used to fix the normalization of the (matter) power 
spectrum. The corresponding as in a given cosmological model then depends on the shape of the 
power spectrum. Analytic approximations to a{R) have already been given by equations (^) and 



(14). Using R = 8h ^ Mpc gives the following approximation for as- 

o_oosa-n)/2 



= 0) 



^8 



0.0136 + 0.294 /i0-803 + 0.109/i2 ' 
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1 + ^5.110.116(0^ 

1 + 0.116(17^ hy 



as{n,) = asin, = 0) , ;n.;..o\Lo 893 ' (15) 



where Qis = Qrms-Ps/18 ^fK. 



The masses and abundances of rich clusters of galaxies provide an independent constraint 
on as- Evrard (1989) obtained erg ~ 0.5 — 0.7 in the standard CDM model from high velocity 
dispersion clusters. Henry & Arnaud (1991) obtained cjg = 0.59 ± 0.02 for scale-free spatially-flat 
models from their X-ray cluster temperature function. White, Efstathiou and Frenk (1993) found 
CTg ~ 0.52 — 0.62 for the standard CDM model using the Henry & Arnaud cluster data and the 
Press-Schechter (1974) approximation (assuming a linear overdensity of 5c = 1.68). Liddle et al. 
(1996) used the same data and also took into account the redshift dependence in the relation 
between cluster virial masses and temperatures. They considered both open and flat CDM models 
and found = 0.59^'^g after allowing 1.6 < Sc < 1.8. 

Here I use the Press-Schechter approximation (1974) to investigate the constraint imposed 
on the neutrino fraction in CDM+HDM models by the cluster data. The comoving number 
density of objects of mass above M at redshift z is given by 



N{M,z) = dlnM' 

J\nM 



2 po{z) 



TT M' a{M',z) 



dlno" 



din M' 



exp 



2cj2(M',z) 



(16) 



where po{z) denotes the background mass density of the universe, and 5c is a free parameter 
characterizing the linear overdensity at the onset of gravitational collapse. The dependence on 
cosmological models is via cr(M), which is given by equations (|l|) and (|l|) with M = A-kR^ po/3. 

White et al. (1993a) approximate a{M) in the CDM model as a single power-law for the 
cluster mass range. As Figure 5 illustrates, this is not a good approximation for CDM+HDM 
models because cj{M) turns over more sharply around lO^^-lO^^M© due to the smaller power in 
P{k). Liddle et al. (1995, 1996) adopted an improved approximation as a function of T = Vlmh for 
open and flat CDM models, but they did not give a formula for o"(r2i^) for CDM+HDM models. To 
compute the number density N{M), one should use either a computed from P{k) in equation (11) 
or the analytic a given by equations (^) and (|l^ above. 

Figure 8 shows the present-day number density of cluster-scale objects as a function of mass 
for various models. Dependence of the cluster number density on the overdensity 6c, normalization 
Qrms-PS) and spectral index n is also illustrated. The sensitivity to the parameters arises because 
cluster masses reside in the Gaussian tails of the Press-Schechter curves. Our experience with 
A^-body simulations is that the best-fit 6c can range from 1.68 to 1.8 (Ma & Bertschinger 1994), 
so both values are shown for comparison. 

Figure 9 illustrates two effects that would help to lower the predicted A^(> M): decreasing 
Qrms-PS or increasing Qf,. The value Qrms-PS = 16.4 pK is the lower 1-a normalization given 
in Bennett et al. (1996). The high value of 0^ = 0.1 (for h = 0.5) is preferred if the recent 
measurement of the deuterium to hydrogen ratio in a quasar absorption line system (Tytler, Fan, 
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k, Buries 1996) is accurate and reflects the primordial value. Another possibility not studied here 
is dividing the contribution of massive neutrinos to Vly among two or more species (Primack et al. 
1995). For a fixed ^7^, the longer free streaming distance of the lighter neutrinos helps to lower the 
power spectrum slightly and therefore gives a lower A^(> M). 

The two data points in Figures 8 and 9 indicate the number densities of clusters with X-ray 
temperatures exceeding HbT = 3.7 and 7 keV from Henry & Arnaud (1991). The mass range 
shown for the ksT > 3.7 keV clusters is taken from White et al. (1993a), where the upper 
limit 5.5 x IO^^/i^-'^Mq was estimated from cluster velocity dispersions assuming an isotropic 
distribution, while the lower limit 4.2 x IQ^^hT^MQ was converted from the X-ray temperature 
assuming an isothermal gas and a density profile of p oc r~^. Both values correspond to masses 
within one AbcU radius 1.5 /i^^ Mpc and are extrapolated from the observed values at ~ 0.5/i~^ 
Mpc. The mass range plotted for the ksT > 7 keV clusters is from the more recent analysis 
of Liddle et al. (1996), who used the hydrodynamical simulations by White et al. (1993b) to 
calibrate the cluster masses within one Abell radius for a given X-ray temperature. Using the 
density profile p oc r~^-^ (White et al. 1993b), Liddle et al. then converted the mass above to 
the virial mass 

{l.2toi) X lO^S/i-^M© for a cluster with a mean X-ray temperature of 7 keV in a 
critical-density universe. 

In comparison with these published values. Figure 8 indicates that all Qrms-PS = 18 fJ-K, n = 1 
models with ^ 0.3 and lit = 0.05 overpredict the cluster number density for 5c = 1.68-1.8. 
Models with n = 0.9 — 0.95 agree better, with the preferred range of dependent on the value 
of Sc and whether or not there is a tensor contribution. Figure 8 shows the same tilted models as 
in Figure 6. For n = 0.9 and no tensor mode (normalized to Qrms-PS = 19.2 /xK), the preferred 
range is 0.2 ^ ^ 0.4, or 0.65 ^ as ^ 0.8, for 6c = 1.68 - 1.8. For n = 0.95 and T/S = 7(1 - n) 
(with Qrms-PS = 18.5 IjlK), the allowed range is 0.1 ;^ $7,/ ;^ 0.3, or 0.6 ^ ag ^ 0.7. This preference 
for n ~ 0.9 — 0.95 is consistent with the conclusion drawn from the reconstructed power spectrum 
in the previous section. A lower Qrms-PS or higher fi^ in n = 1 models also helps to decrease the 
cluster number density, but the effects are smaller (see Fig. 9). 

It should be kept in mind that the mass ranges adopted by White et al. (1993a) and Liddle 
et al. (1995, 1996) are obtained under specific assumptions. Cluster mass determination is by 
no means a settled issue. Temperature gradients, substructure, and ellipticity can all affect the 
X-ray mass estimates (see, e.g., Tsai, Katz, Sz Bertschinger 1994). Cluster mass function has 
also been estimated from velocity dispersions of cluster galaxies (e.g., Biviano et al. 1993), but 
inclusion of substructure in the analysis lowers the normalization (Bird 1995). We can invert the 
problem and ask what cluster mass is needed for the currently favored n = 1 CDM+HDM models 
to produce the observed cluster abundance. Within the narrow window Cli, ~ 0.1 - 0.2 allowed 
by both the low- and high-redshift constraints (Ma & Bertschinger 1994; Ma 1995), Figure 8 
indicates that the 3.7 and 7 keV clusters would have to have a median mass of ~ 1 - 1.5 X 10^^ 
and ~ 2 — 4 X 10^^h~^ Mq, respectively, for Q^, = 0.1 - 0.2 models (n = 1) to produce the 
observed cluster abundance. This seems very unlikely, given that the hottest and most luminous 
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X-ray cluster Abell 2163 {ksT = 12 — 15 keV) is recently estimated to have a mass of only 
Ma = 1.07 ±0.13 X W^^H-^Mq out to the Abell radius (Markevitch et al. 1996). For a power-law 
density profile p(r) oc r^" at the outer part of the cluster, this translates to a virial mass of 
My PS 1.5(3/°-^)Myi, that is, only about 1.2 - 1.35 x IO^^H-'^Mq for a ranging from 2.4 from 
simulations (White et al. 1993b) to the isothermal a = 2. 

4. Conclusions 

This paper presents numerical results and analytical approximations to the linear power 
spectrum and the density field growth rate for a class of CDM+HDM models with neutrino 
fraction ^ 0.3 (neutrino mass rrii, ^ 7 eV) that are currently of much cosmological interest. 
The models with $7^ = (i.e. the standard CDM) and ^ 0.3 have difficulty explaining various 
observations, leaving the range of the neutrino masses studied in the paper as the only hope for 
the standard CDM+HDM models. The analytical functions given by equations (|^)-(|ll|) and ( |l3[ ) 
and (U) in this paper give accurate approximations to the linear growth rate, the separate cold 
and hot power spectra, the density-weighted power spectrum, and the rms mass fluctuations for 
k ^20h Mpc-\ z ^ 15, and 0.5 ^ h ^ 0.8 in CDM+HDM models with ^ 0.3. They should 
be useful as the input for most linear calculations and for the initial conditions of numerical 
simulations of structure formation in these models. 

It is shown that the difference in the growth rate and the power spectrum between the 
CDM+HDM and the standard CDM models can be entirely characterized by the free-streaming 
distance of the neutrinos. A shape parameter T^ = a^^'^Qi, h?' is introduced to explain the 
behavior of the power spectra in Figure 1; it also helps to greatly simplify the functional forms of 
equations (^), (P), and ([To|), and to improve the accuracy of the fits. 

Linear calculations with these power spectra as input are performed to study the dependence 
of cluster abundance on the spectral index n, normalization Qrms-PSi tensor-to-scalar ratio T/S, 
baryon fraction $7^, and neutrino masses in Vti, ^ 0.3 models. The theoretical predictions are 
compared to the reconstructed P{k) by Peacock &; Dodds (1994) and the X-ray cluster abundance 
and masses from Henry &; Arnaud (1991). These tests are chosen because the linear theory (with 
proper calibration for the clusters) should be adequate on such large scales. When normalized 
to the 4-year COBE Qrms-PS = 18/iii', all n = 1, fi;, = 0.05, VL^, <^ 0.3 models predict too much 
power by up to a factor of ~ 2 when compared to Peacock & Dodds (1994) and Herny & Arnaud 
(1991). A lower Qrms-PS or higher Ofe both helps to reduce the discrepancy. Also considered is a 
slight tilt of n = 0.9 — 0.95 in the power spectrum, which brings the theoretical predictions into 
good agreement with data. 

Due to the heterogeneous nature of the Peacock & Dodds sample and the uncertainties in 
cluster abundance and mass, it is difficult to attach a meaningful statistical significance to the 
validity of the various models studied in this paper. It should also be kept in mind that the 
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conclusion that the n = 0.9 — 0.95 models provide a closer match to the data depends on several 

implicit assumptions. These include the scale-independent bias used in Peacock & Dodds (cf. 
Cen &: Ostriker 1992) and the isothermal and spherical model for clusters used in White et al. 
(1993a) and Liddle et al. (1996) to convert X-ray temperature into virialized mass. These are 
obviously simplifications, but modeling the deviations from such assumptions would require the 
introduction of parameters that can not be well determined by current data. Nonetheless, the 
theoretical calculations carried out in this paper should provide useful predictions that can be 
tested by future data. 

The author is grateful to Ed Bertschinger and Peter Goldreich for valuable discussions, Una 
Hwang, Mark Metzger, HouJun Mo, Paul Steinhardt, and Ned Wright for helpful comments, and 
K. Gorski and Richard Mushotzky for useful information. Parts of the Boltzmann calculations 
were performed at the National Center for Supercomputing Applications. Support from a PMA 
Division Fellowship at Caltech is acknowledged. 
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Table 1. Summary of Analytic Approximations 



Symbol 


Equation 


Coefficients 1 


2 


3 


4 


5 


/ 


dD 


0.1161 


1.363 








Pci^u = 0) 




2.34 


3.89 


16.1 


5.46 


6.71 






2.205 


4.05 


18.3 


8.725 


8.0 




® 


0.01647 


2.803e-5 


10.9 


3.259 






m 


0.0015 


-0.1207 


0.1015 


-0.01618 


0.001711 


P{^u) 




0.004321 


2.217e-6 


11.63 


3.317 




= 0) 


m 


0.01359 


0.05541 


0.001702 


0.8032 






(0) 


0.7396 


-0.8927 


5.106 







^For zero baryons, or modify variable T in eqn. (|^) for general baryon fraction (see 
text) 

'^High accuracy fit (with ^ 1% error) for the standard CDM model with 5% baryons 
■^Density-weighted spectrum P = {Q^P^^"^ + (1 - ^i^)Pc^'^}'^ 
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Fig. 1. — Present-day linear power spectra for the standard CDM model (dotted) and four 
CDM+HDM models with n = 1, Qb, and Hq = 50 km Mpc~^, computed from integrations of 
the Boltzmann equations. The four CDM+HDM models have = 0.05,0.1,0.2, and 0.3, and the 
power in the cold (solid) and hot (dashed) components are shown separately. All are normalized 
to the 4-year COBE result Qrms-PS = 18 jiK (Gorski et al. 1996). 
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Fig. 2. — Direct integration of the Boltzmann equations (solid) vs. analytic fitting results (dashed) 
for the present-day CDM and HDM power spectra and the growth rate in four CDM+HDM models 
(n = 1,Hq = 50) with Jljy = 0.05,0.1,0.2, and 0.3. (a) Growth rate f{k,a,il.^) of the density field. 
The fitting formula is given by eq. (P). (b) Ratio of the CDM power spectrum in CDM+HDM 
models to that in the pure CDM model, g = Pc{0,u) / Pc{0,i, = 0). The fitting formula is given by 
eq. (H). (c) Ratio of the HDM to CDM power spectrum, H = P^in^) / Pc{n^) . The fitting formula 
is given by eq. (p^. 
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Fig. 3. — Same as Figure |^ but for a = 0.1 
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Fig. 4. — Scaling of k with the Hubble constant for functions /, g and TC in the = 0.2, n = 1 
CDM+HDM model at a = 1. As discussed in the text, / and H scale perfectly with k/h^, while 
k/h? IS a, good approximation for g for a large range of k. 
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Fig. 5. — Root-mean-square of the linear mass fluctuations a{R,^i^) in spheres of radius Rh~^ 
Mpc for n = 1 and Hq = 50 models with = 0.0,0.05,0.1,0.2, and 0.3 (top down). The solid 
and dashed curves are from numerical integration and the fitting formulas eqs. ( [l^ ) and (14), 
respectively. The dotted curves show the scaled a~^a{R, , a) for a = 0.1, which start to deviate 
from a{R,Vty,a = 1) only at small scale. All models are normalized to Qrms-PS = 18 at large 
R. 
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Fig. 6. — Density-weighted linear power spectrum for four models with neutrino fractions 
= 0.0,0.1,0.2 and 0.3 (solid curves from top down at large k), compared with the linear 
power spectrum reconstructed by Peacock & Dodds (1994) from galaxy and cluster surveys (filled 
symbols). The baryon fraction is O;, = 0.05 and Hq = 50 km Mpc~^. The bottom two panels 
show the excellent agreement of CDM+HDM models with the data when the spectral index is 
reduced to n = 0.95 with tensor fluctuations (middle) or to ra = 0.9 without tensor fluctuations 
(bottom). The 4-year COBE result is used for the normalization, and the dependence of ^5rms— PS 
on n is taken into account (see text). The dotted curves show an n = 1 CDM-type spectrum with 
a shape parameter of F = Qh = 0.25. 
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Fig. 7. — Same as Figure ^ but for Hq = 65 km Mpc~^ (with O;, = 0.05). A larger tilt of 
n ~ 0.9 (with tensor) or n ~ 0.8 (without tensor) is needed to reduce the excess power in the n = 1 
model. 
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Fig. 8. — Present-day comoving number density of cluster-mass objects as a function of mass for 
the standard CDM and three CDM-I-HDM models. Different panels illustrate the effect of varying 
the spectral index n and the parameter 5c in the Press-Schechter approximation. Top panels are for 
n = 1 models with Qi, = 0, 0.1, 0.2, and 0.3 (from top down), normalized to the 4- year COBE result 
Qrms-PS = 18 fiK (Gorski et al. 1996). The solid and dotted curves in the upper left panel compare 
A^(> M) computed from numerical vs. fitted a. Bottom panels show the same tilted models with 
= 0, 0.1, 0.2, and 0.3 (from top down) as in Figs. ^ and |^ The dependence of Qrms-PS on n is 
taken into account: Qrms-PS = 19.2 fiK for n = 0.9 and Qrms-PS = 18.5 for n = 0.95. The data 
points indicate the number density of clusters (Henry & Arnaud 1991) with X-ray temperature 
exceeding 3.7 (upper left) and 7 keV (lower right), respectively. The corresponding mass range is 
taken from White et al. (1993a) and Liddle et al. (1996). 
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Fig. 9. — Effects of changing Qrms-PS (upper panels) and Qb (lower panels) on the number density 
of cluster-mass objects. For clarity, only fi^ = and 0.3 models are shown. All models have n = 1 
and h = 0.5. The top panels assume Qi, = 0.05; the bottom panels assume Qrms-PS = 18 jJ-K. The 
data points are the same as in Fig. & 



